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ABSTRACT. 

This paper presents implementation of new algorithm Particle Swarm Optimization (PSO)for Energy 
Saving through minimizing power losses. The PSO Algorithm Solution is tested in standard IEEE 30 
Bus system. The objective is to optimize the reactive power dispatch with optimal setting of control 
variables without violating inequality constraints and satisfying equality constraint. Control Variables 
are of both types: Continuous and Discrete. The continuous control variables are generator bus voltage 
magnitudes; whereas the discrete variables are transformer tap settings and reactive power of shunt 
compensators (Capacitor banks) . 
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I. INTRODUCTION 

Power system economics procedure includes of two aspects: active power regulation and reactive 
power dispatch. This forms a global optimization problem of a large-scale industrial system.The reactive power 
problem is less manageable to solve than the active power problem due to its more complicated relationship 
between variables. The reactive power problem is largelyassociated to voltage stability. Reactive power and 
voltage control is incredibly essential for the right operation and control of power system. Reactive power 
dispatch is one of the necessary tasks in the operation and control of power system. Voltage stability is a 
drawback in power systems that are heavily loaded, faulted or have a deficiency of reactive power. The 
character of voltage stability may be analysed by examining the generation, transmission and consumption of 
reactive power. Transfer of reactive Power is toughbecause of extremely high reactive power losses; that's why 
the reactive power needed for voltage control is generated and consumed at the control area. Optimization is a 
mathematical procedurewhich discusses the finding of maxima or minima of functions in some realistic region. 
There's no industry or businessthat is not involved in solving optimization problems. By Optimizing reactive 
power Dispatch in Power systems, the maximum active power transfer capability to the distribution systems can 
be improved. Stand-by reactive power sources (capacitor banks generally) are needed for loss minimization, in 
order to maintain the voltage stability in the Power systems [11-12]. 

For solving all optimization problems, there is no known single optimization method available. For 
solving the differentkinds of optimization problems, plenty of optimization techniques have been established in 
recent years. Linear programming (LP), non-linear programming and gradient based techniques were traditional 
optimization techniques [1], [4] for solving Reactive Power optimization problems. Since, Approximations are 
used in linearized models, thus LP results don't signify optimal result for objective function utilized in reactive 
power optimization problem. Traditional solution strategies have tendency to converge to a local optimal 
solution instead of the global one. Expert System methodologies [5] have been recommended for the reactive 
power based calculations. Expert System methodology is based mostly on 'if-then' based rules. Evolutionary 
computational techniques like Genetic algorithm (GA), Evolutionary programming (EP) and Evolutionary 
strategy have additionally been projected to solve the optimizations problems relating to the reactive power 
dispatch [6-8]. The contemporary (non-traditional) optimization approaches are very powerful and popular 
approaches for solving complex engineering problems. These approaches are neural networks, genetic 
algorithm, ant colony optimization, fuzzy optimization and particle swarm optimization algorithm Particle 
Swarm Optimization (PSO) stands as acomparatively new, modern, and powerful technique of optimization that 
has been virtually shown to perform well on several of these optimization problems [13-14]. PSO exists as a 
population based stochastic optimization technique. PSO algorithm is applied while not violating inequality 
constraints and satisfying equality constraint.The aim of minimising reactive power losses is achieved by 
appropriate adjustment of reactive power variables like generator voltage magnitudes (V gi ), reactive power 
generation of capacitor banks (Q ci ) and transformer tap settings (t k ) [7-10]. In electrical power system, Reactive 
Power Loss Minimization problem is taken into account as a static, non-linear, single objective optimization 
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problem. The suggested PSO algorithm solution has been experimented on the standard IEEE 30-Bus test 
system with both continuous and discrete control variables despite the fact that keeping the system under safe 
voltage stability limit. The recommended algorithm shows better results 

II. PROBLEM FORMULATION 

The customary optimization problem is often written in the following form, 
Minimise F(x) (the objective function) 
subject to: 

hi(x) = 0, i = 1, 2, n (equality constraints) 
g/x) = 0, j = 1, 2, ...,m (inequality constraints) 

The reactive power optimization problem targets to minimize the power losses in the transmission 
network and improve voltage profile while satisfying the unit and system constraints. The aim is achieved by 
appropriate setting of reactive power variables like generator voltage magnitudes (V gi ), reactive power 
generation of capacitor banks (Q ci ) and transformer tap settings (t k ) [7-8]. The equality constraints are 
power/reactive power equalities, the inequality constraints consist of bus voltage constraints, generator reactive 
power constraints, reactive power capacity constraints and the transformer tap position constraints, etc. At this 
juncture the reactive power dispatch problem is treated as a single objective optimization problem by linear 
combination of two objective functions i.e. I^ 33 

F = mm.P iQSS [1] 

A. Energy Saving through Minimization of power system losses (Pi oss ) 

The RPD problem targets at saving of energy by minimizing the real power loss in a power system while 
satisfying the unit and system constraints. This objective is accomplished by appropriate adjustment of reactive 
power variables like generator voltage magnitudes (V G i), reactive power generation of capacitor banks (Qa) and 
transformer tap settings (T k ). 

The minimization of system real power losses (MW) is calculated as follows: 

□l 

Min F = H„ fl = ^ gJVi* + If - ZViVj cosfe - a,)] 

The real power loss given by (Pl oss ) is a non-linear function of bus voltages and phase angles that are a function 
of control variables and ni is the number of transmission lines; g k is the conductance of the k th line; Vi and Vjare 
the voltage magnitude at the end buses i and j of the k th line, respectively, and £ £ and are the voltage phase 

angles at the end buses i and j. 

III. CONSTRAINTS 

The real power loss (given by equation) is treated as a non-linear function of bus voltages and phase 
angles that are functions of control variables. The minimization problem is subjected to the following equality 
and inequality constraints: 

Equality constraints 

These constraints are typical load flow equations which can be represented as follows 

1. Real Power Constraints: 

Pa ~ Cfai " ViSj^Vj (tyCosfe - £j) + B^Sinfe - 6,)) = 0 ...[2] 

I = 1,2 Ng 

2. Reactive Power Constraints: 

ttci " Qoi -ViSj^Vj^anfe - B^Cosfe = 0,...[3] 

t = ^PQ 

Where, 

Vi= Voltage magnitude at bus I 
Vj = Voltage magnitude at bus j 

Qj = Real and reactive powers injected into network at bus i 
Gjj, Ejj = Mutual conductance and susceptance between bus i and bus j 
Q^i = Reactive power generation at bus i 
N B — 1 = Total number of buses excluding slack bus 
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iV^ t =Number of PQ buses 

6ij = Voltage angle difference between bus i and bus j 

Inequality constraints : 

3. Bus Voltage magnitude constraints: 
V-™ 11 < Vi < V- maj: ; ibN b ... [4] 

4. Transformer Tap position constraints: 
tj^<t k <t™* ; «t T ...[5] 

5. Generator bus reactive power constraint: 
Q^T<^ < $%T' r id* ff ..[6] 



6. Reactive power source capacity constraints: 

7. Transmission line flow constraints: 

ki <ur*i....[8] 
uj <\sr*\...[9] 

8. Generation capacity constraint: 
Pg* < Pci < Eg" ; iEN B [10] 

The total power generation should cover the overall demand P D and the real power loss in transmission lines 
P L . This relation is often expressed by Power Balance Constraint: 

9. Power balance constraint: 

Sffl^L =P B +P L ....[U] 

The symbols used are as follows: 
t]<= Tap setting of transformer at branch k 

Qci= Reactive power generated by 1 th capacitor bank 

Q^i = Reactive power generation at bus i 

Si = Apparent power flow through the I th branch 

N B = Total number of buses 

g k = Conductance of buses 

N T = Number of tap-setting transformer branches 

N c = Number of capacitor banks 
N g = Number of generator buses 

The control variables for voltage-control problem, which will be modified by the Particle Swarm 
optimization process, are: 
a. Voltages magnitude at voltage-controlled buses (PV-buses) including the slack bus. 

b. Transformers tap settings. 

c. Adjustable shunt capacitor banks. 

IV. PARTICLE SWARM OPTIMIZATION 

PSO is an acronym for Particle Swarm Optimization. Particle Swarm Algorithm was introduced by 
Kennedy and Eberhart in 1995 [14]. PSO is a swarm intelligence method for global optimization.Particle Swarm 
Optimization is a concept introduced for the optimization of nonlinear functions using particle swarm 
methodology. Basically Particle Swarm Optimization is a method for optimization of continuous nonlinear 
functions. The method was discovered through simulation of a simplified social model. Particle Swarm 
Optimization comprises a very simple concept, and paradigms can be implemented in a few lines of computer 
code. It requires only primitive mathematical operators, and is computationally inexpensive in terms of both 
memory requirements and speed. Early testing has found the implementation to be effective with several kinds 
of problems. 
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PSO is based on the natural process of group communication to share individual knowledge when a 
group of birds or insects search food or migrate and so forth in a searching space, although all birds or insects do 
not know where the best position is. But from the nature of the social behaviour, if any member can find out a 
desirable path to go, the rest of the members will follow quickly. PSO traces its evolution to the emergent motion 
of a flock of birds searching for food. PSO uses a number of particles that constitute a swarm. Each particle 
traverses the search space looking for the global minimum (or maximum). In a PSO system, particles fly around 
in a multidimensional search space. During flight, each particle adjusts its position according to its own 
experience, and the experience of neighbouring particles, making use of the best position encountered by itself 
and its neighbours. The swarm direction of a particle is defined by the set of particles neighbouring to the 
particle and its history experience. 

The basic elements of the PSO techniques are defined as: 

1. Particle X (t) : It is a candidate solution described by a k- dimensional real-valued vector, where k is the 
number of optimized parameters. At iteration i, the jth particle X (i,j) can be expressed as: 

Xi(t)=[xi,i(t); Xi, 2 (t); ;xi, k (t)]. 

Where: x's are the optimized parameters and d signifies number of control variables 

2. Population: It is basically a set of n particles at iteration i. 
Pop (i) = [ Xi®, X 2 (i), X n (i)] T 

Where n signifies the number of candidate solutions 

3. Swarm: Swarm is defined as an apparently unsystematic population of moving particles that tend to bunch 
together while each particle appears to be moving in a random direction. 

4. Particle velocity V (t):Particle velocity is the velocity of the moving particles signified by a d-dimensional 
real valued vector. Particle Velocity is the step size of the swarm. It is the velocity of the moving particles 
represented by a k-dimensional real- valued vector. At time t, the i th particle Vi(t) can be described as 

Vi (t)=[Vi,i(t); Vi, 2 (t); ;vi, k (t)]. 

5. Inertia weight «(t): It is a regulation parameter, which is used to regulate the impact of the past (previous) 
velocity on the present velocity. Hence, it effects the trade-off between the global and local exploration 
capacities of the particles. For the initial stages of the search method, large inertia weight to reinforce the global 
exploration is usually recommended while it must be reduced at the last stages for higher local exploration. 
Therefore, the inertia factor drops linearly from about 0.9 to 0.4 throughout a run. In general, the inertia weight 
factor is set according to the equation given below: 

W = : x iter 

All the control variables transformer tap positions and switch -able shunt capacitor banks are integer variables 
and not continuous variables. Therefore, the value of the inertia weight is considered to be 1 in this study. 

6. Individual best X* (t): When particles are moving through the search space , it matches its fitness value at 
the existing position to the best fitness value it has ever grasped at any iteration up to the current iteration. The 
best position that is related with the best fitness faced so far is called the individual best X* (i). For every 
particle in the swarm, X* (i) may be determined and updated throughout the search. 

7. Global best X** (t): It is the best position among all of the individual best positions achieved so far. 
Various steps concerned with the implementation of PSO to the RPO problem are: 

Step 1: Firstly scan the Input parameters of the system (bus, line and generator data) and alsoidentify the lower 
and upper boundaries of every variable. For N generators, optimization is applied out for N-l generators and 
generator of maximum capacity is considered at slack bus. 

Step 2: Then the particles of the population are randomly initialized i.e. are randomly selected between the 
respective minimum and maximum values. Also assign the velocity V initially between [-1 and 1]. 
Step 3: Obtain power flow solution and compute losses by Newton -Raphson method. 
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Update Counter 



Compute the 
Objective Function 



i 



Update Velocity 



I 



Update Individual best 
fPbid if weeded 



Update Position 



I 



Update Global best 
(g*md if needed 



Xo 




Optimal Control Variable 
Setting 



Fig.: Flowchart for Basic PSO Algorithm 

Step 4: The best fitness is assigned aspj eB t. At this stage thepj ei5 tis also theg^ est . 
Step 5: Iteration^ i+1 is updated. 

Step 6: Modify the inertia weight w given by W = ^ lVmiLE x iter 



Lfsr, 



■max 



Step 7: Update the velocity v of each particle according to the stated equation 

V<JkJ f i+l} = rPtt^Ot Crrand^ptexAhk) - JTCftJ.O) +■ C2 Vanri * (fe) - JT(fc,;, 0) 

Step 5: Position of each particle is also modified as per the mentioned equation. If a particle violates the position 
limits in any dimension, its position is set at the right limit. 
x(k f j f il) x(kj f l0 v{k..j., 0 
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Table: optimal parameter setting for PSO 



Parameters 


Number of iterations 


300 


Cognitive constant, ci 


2.0 


Social constant, C2 


2.0 


Max. and Min. inertia weights W 


0.4 and 0.95 


Population size 


30 



Step 9: Evaluation of each particle is done according to its updated position by runningpower flow and calculate 
the fitness function. If the evaluation value of each particle is better than the previous Pb &s1 then the current 

value is set to be Pse-st • If the best Pb bs1 is better than gs&st » the value is set to be gs&st- 

Step 10: If one of the stopping criteria is fulfilled then we go to Step 11. Else, we go to Step 5. 
Stepll: g Best is the optimal/best value that is newest generated by the particle. 




Fig.: IEEE 30 Bus System 

V. RESULTS 
Energy Saving through Minimization of systempower losses (Pi 0S s) 

The proposed algorithm is run with minimization of real power losses as the objective function. As mentioned 
above, the real power settings of the generators are taken from [15-16]. The algorithm reaches a minimum loss 
of 5.3191MW. IEEE30 bus system is shown. The optimal values of the control variables are given in table 
shown above. . 



Table 1: Test results of proposed approach. 



Proposed Method 


Power Loss (MW) 


Particle Swarm 
Optimization 


5.3191 





Min 

• 


Max. 


Initial(Ba 

se case) 


Proposed 
PSO 

algorithm 


VI 


1.0 


1.1 


1.05 


1.0824 


V2 


1.0 


1.1 


1.04 


1.0470 


V5 


1.0 


1.1 


1.01 


1.0347 


V8 


1.0 


1.1 


1.01 


1.0209 
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Vll 


1.0 


1.1 


1.05 


1.0376 


V13 


1.0 


1.1 


1.05 


1.0402 


Til 


1.0 


1.1 


1.078 


1.0196 


T12 


1.0 


1.1 


1.069 


1.0783 


T15 


1.0 


1.1 


1.032 


1.0573 


T36 


1.0 


1.1 


1.068 


1.0963 


QclO 


0.0 


5.0 


0.0 


1.2677 


Qcl2 


0.0 


5.0 


0.0 


1.0610 


Qcl5 


0.0 


5.0 


0.0 


0.8607 


Qcl7 


0.0 


5.0 


0.0 


0 


Qc20 


0.0 


5.0 


0.0 


2.5792 


Qc21 


0.0 


5.0 


0.0 


1.7678 


Qc23 


0.0 


5.0 


0.0 


1.6902 


Qc24 


0.0 


5.0 


0.0 


0.5076 


Qc29 


0.0 


5.0 


0.0 


0.6881 


Power loss(MW) 


5.8708 


5.3191 



Table 2:- Best results of individually run of pi oss as main function (IEEE- 30 Bus) 

Total Energy Saving : 

From the table (shown below), 

Reduction in losses - (5.8708 - 5.3191) = 0.5517 MW = 551.7 KW 

Converting these reduced active power losses in form of energy we find: 

Saved Energy in One Hour => 551.7 KWh 

Saved Energy in One Day => 551.7 *24 = 13240.8 KWh 

Saved Energy in One Week => 551.7 *24*7 = 92685.6 KWh 

Saved Energy in One Year - 551.7 *24*365 = 4832892 KWh 

So for an IEEE 30 bus system, estimation of energy saving is shown. For a standard IEEE 30 Bus system, 
4832892 KWh (Units) is saved in a year by using PSO. 

VI. CONCLUSION 

A new improved integer coding Particle Swarm Algorithm is presented to solve this problem. The main 
objective is to minimize the active power loss in the network, while satisfying all the power system operation 
constraints. The particle swarm algorithm has been coded as well as the power flow fast-decoupled method 
using MATLAB. The simulation results show that PSO algorithm always leads to a better result. 
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